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Abstract 

For an infinite system of nucleons interacting through a central spin-isospin schematic 
force we discuss how the Hartree-Fock theory at finite temperature T yields back, 
in the T = limit, the standard zero-temperature Feynman theory when there is no 
symmetry breaking. The attention is focused on the mechanism of cancellation of 
the higher order Hartree-Fock diagrams and on the dependence of this cancellation 
upon the range of the interaction. When a symmetry breaking takes place it turns 
out that more iterations are required to reach the self-consistent Hartree-Fock so- 
lution, because the cancellation of the Hartree-Fock diagrams of order higher than 
one no longer occurs. We explore in particular the case of an explicit symmetry 
breaking induced by a constant, uniform magnetic field B acting on a system of 
neutrons. Here we compare calculations performed using either the single-particle 
Matsubara propagator or the zero-temperature polarization propagator, discussing 
under which perturbative scheme they lead to identical results (if B is not too large) . 
We finally address the issue of the spontaneous symmetry breaking for a system of 
neutrons using the technique of the anomalous propagator: in this framework we 
recover the Stoner equation and the critical values of the interaction corresponding 
to a transition to a ferromagnetic phase. 

Key words: Nuclear matter, finite temperature, Hartree-Fock, random phase 
approximation 
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1 Introduction 



In this work we first discuss the problem of the Hartree-Fock (HP) mean field 
for an infinite system of fermions (nucleons) at finite temperature T in the 
context of Matsubara theory. In particular we study under which conditions 
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and how the results obtained in the Matsubara framework in the T — > hmit 
coincide with those of the T — Feynman theory. 

As is well-known, the two formalisms lead to identical results in the absence 
of symmetry breaking, according to the old theorem for spin one-half fermions 
of Kohn, Luttingcr and Ward (KLW) [1,2], at least as far as the ground state 
energy is concerned. Here we extend the theorem to the singlc-particlc prop- 
agator (and hence to any one-body observable), aiming at transparently dis- 
playing how the cancellation of the diagrams, the mechanism at the basis of 
KLW, is reahzed and how the exact T — limit is retrieved. 

Specifically we shall assess quantitatively the magnitude of the diagrams that, 
while canceling in the T — > limit, actually substantially contribute at finite 
T (not, however, in finite systems [1,2,5]). Concerning the size of their con- 
tribution to observables like the chemical potential and the magnetization, it 
turns out to be larger in the proximity of the Fermi temperature Tp owing to 
the temperature dependence of the self-energy. The latter, in fact, displays a 
marked change near Tp reflecting the transition from a quantal to a classical 
regime. This transition is sensitively affected by the range of the interaction: 
hence in our study we employ both a finite and a zero-range force, however 
of schematic nature since we have no pretense of performing a realistic cal- 
culation, but rather we aim to investigate the generic temperature behavior 
of observables that are relevant for the finite temperature physics (such as 
the chemical potential and the magnetization) in the absence or presence of a 
symmetry breaking. 

We address this last issue in the second part of the paper. When a symmetry 
is broken then KLW no longer holds and accordingly we explore the impact 
on the HF field of such an occurrence, both when the breaking is induced 
explicitly by an external field (as an example we shall consider the case of a 
constant, uniform magnetic field acting on a system of neutrons) and when 
it is spontaneous. In the first instance we show how the diagrams that would 
cancel each other at small temperature in a situation of pure symmetry no 
longer do so: actually their contribution grows, at T — > 0, keeps the value 
attained a,t T — Tp if the applied magnetic field is large enough. 

Moreover, and interestingly, in this situation the successive iterations approach 
the self-consistent mean field solution smoothly or oscillating around the HF 
value depending upon the sign of the two-body interaction among neutrons. 
A ferromagnetic force gently leads to the HF mean field, whereas an antifer- 
romagnetic force entails an approach to the latter oscillating around its value, 
the magnitude of the oscillations increasing with the strength of the interac- 
tion. This behavior contrasts the one occurring in the absence of symmetry 
breaking where the HF solution is always smoothly reached and, in fact, it is 
also rapidly reached — at least for a pure exchange interaction (the one we 
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confine ourselves to consider) — at small and large temperature: in the first 
case because of the KLW theorem, in the second because the large T domain 
is where quantum physics is gone and classical physics sets in. 

The unusual property, referred to above, of the finite temperature HF solutions 
in presence of a magnetic field B stems in part from the remarkable occurrence 
that B, in breaking the spin rotational symmetry of the system, induces in the 
interaction matrix element a direct term which otherwise would be absent, the 
force we employ being of pure exchange character. This is of importance when 
one aims to recover the vanishing temperature limit of the observables, say the 
magnetization, computed in the Matsubara HF formalism in the framework 
of the T = Feynman theory. This in fact turns out to be possible, in spite 
of the breakdown of the KLW theorem, using linear response theory, whose 
applicability is however limited to a specific range of B, which we are able to 
assess quantitatively. 

Concerning instead the case of spontaneous symmetry breaking, the non- 
relativistic quantum field theory describes its occurrence at T = through 
the linearization of the equations of motion, i.e. at the mean field level, if the 
interaction is strong enough [3] . This approach coincides with the many-body 
Hartree mean field theory at T = 0, which indeed exhibits a spontaneous 
breaking of symmetry — again for a force strong enough — provided the vari- 
ational search for the determinant yielding the minimum energy is allowed to 
span all of the possible symmetry configurations (we refer to this approach as 
the generalized T — HF) . This is the path we follow in the last Section in 
the framework of the anomalous propagator technique [4] at T = 0. 

Likewise, also the temperature HF theory displays a spontaneous symmetry 
breaking. We prefer, however, to search for the onset of the latter in the frame- 
work of the T — linear response theory which, in parallel to the temperature 
HF, predicts, either in ring or in the random phase approximation (RPA), 
a divergent magnetic susceptibility of the Fermi system when the strength 
of the nucleon-nucleon interaction, assumed to be ferromagnetic, reaches a 
critical value Vent, thus signaling the occurrence of a phase transition. For 
a contact interaction, when the response of the system is computed in ring 
approximation, the value of V^rit thus found is identical to the one obtained 
in [3] or in the generalized Hartree approximation, whereas the RPA response 
yields a lower value (by a factor 2/3), which coincides with the one predicted 
by the T = generalized HF theory. As mentioned above, the same result is 
obtained using the thermal single-particle HF propagator, computed with a 
contact interaction, in the T" — > limit. 

When the symmetry is spontaneously broken, one would like to describe the 
associated Goldstone modes. In connection with the collective excitations of 
the system, at zero temperature it turns out that, for an anti-ferromagnetic 
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(a) (b) (c) (d) 

Fig. 1. The perturbative expansion of the HF irreducible self-energy for a pure 
exchange interaction. 

coupling among the constituents, the Fermi system supports the existence 
of zero-sound collective modes characterized by a linear relationship between 
energy and momentum. These modes correspond to the spin waves (also re- 
ferred to as magnons) found in the antiferromagnetic materials and lower the 
susceptibility. They represent the Goldstone modes, because our system, al- 
though without any cristal structure, still lives in a condensed phase given by 
the antiferromagnetic order of the spin constituents: we accordingly view the 
Fermi gas as paramagnetic at finite T and antiferromagnetic at T = 0. 

On the other hand a ferromagnetic interaction only supports the existence of 
modes embedded in the particle-hole continuum for V < Krit- In this regime 
such a force softens the response function of the system until, when V = Krit, 
the latter diverges at zero excitation energy as the transferred momentum be- 
comes very small. Thus the phase transition is also displayed by the behavior 
of the response function, which basically corresponds to the spin-spin corre- 
lation function. On the other hand, for a ferromagnetic coupling V > Vait 
new collective modes appear, in the direction orthogonal to the one set by 
the spontaneous magnetization of the system. These also represent Goldstone 
modes and will be addressed in a forthcoming paper. 



2 Hcirtree-Fock theory at zero and finite temperature 

In this Section we discuss the mean-field HF theory of a system of nucleons at 
finite temperature. For simplicity wc confine ourselves to consider an infinite 
homogeneous system and a schematic static exchange interaction in the non- 
relativistic limit. 

In this situation it is well-known that at zero temperature the HF mean-field 
theory reduces to first order perturbation theory. Indeed for the single-particle 
wave functions self-consistency is immediately achieved in the first iteration of 
the HF equations owing to the translational invariance of the system. In con- 
formity, the diagrams contributing to the HF self-energy (displayed in Fig. 1) 
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vanish, except for the first one: indeed all of them display in the energy variable 
poles of order > 2 with zero residue. 

For example, considering the diagram (b) of Fig. 1 one has 
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V being the interaction and G° the zero-order propagator in a Fermi gas with 
Fermi momentum kp- Now, in the frequency integral. 
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only the double pole at uji — ui^^ + ir] = hk\/2m + ir] should be considered, 
the contour of the integration lying in the complex upper plane Im a;i > 0. 
Hence one gets 



I=-r, e*K+^'')''^(A;^ _ ^ q for r/ ^ 
and hkewise for all the other diagrams of Fig. 1 of order > 2. 



(3) 



On the other hand, at finite temperature the Matsubara diagrams of Fig. 1 no 
longer vanish. In fact, taking again as an example the second order self-energy, 
one has 
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XV{k-k^)V{k - k2)[G\k,UJn,)fG\k2,UJn,), (4) 



where 



G\k,^n)^- -^0^ — (5) 

is the thermal free propagator at finite T — l/ksP {ks being the Boltzmann's 
constant and ^ the free chemical potential). Evaluating the frequency sums 
according to the standard rule 

hm ^ e--''^°(4, ^nj = lim ^ ''^'^ = (6) 
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and 



e *=! + 1 



2 ' 



(7) 



where c<j„ = {2n + 1)7: / (3h and the sums run over all the positive and negative 
integers, one gets 



h J {2tiY J (27r 



/3(4°'-M) 



e '^i 

^ (0) , ■ (8) 
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The self-energy (8) does not vanish: as a consequence, at finite T, even for an 
infinite system, the HF equations are not trivial. 

We briefly illustrate this point using, as an example, a simple static spin- 
isospin dependent nucleon-nucleon central interaction of the type 

V = di- a2Ti ■ T2v{Xi - X2) (9) 

which, of course, has no pretense of being realistic. 

Since the system is homogeneous, the eigenfunctions are still plane waves, 
whereas the HF eigenvalues read 

^k = 4^ + ^KF{k,n (10) 

where 

n^Uk, T)^-9j j^,npv{k - k') (11) 
is the HF irreducible self energy for the interaction (9), 

nr^ \ (12) 

— * 

the Fermi distribution and v{}i) the Fourier transform of v{x). 

Thus the HF equations at flnite temperature, unlike the case at T = 0, repre- 
sent a non trivial self-consistency problem, because the single-particle energies 
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(the eigenvalues) enter not only in the equations, but in the Fermi distribution 
as well. 



Accordingly, the system of equations 



ej; = ef+h^'kF{k,T) (13a) 



dk' 



(27r) 



3^fc" (13b) 



where the density p of the system is assumed to be fixed, should be solved. 
We have accomplished this task through numerical iterations of the above 
equations (the unknowns being eg and n) until reaching self-consistency (see 
Subsection 2.3 for the numerical results). 

At zero order eg = el^-* and one fixes = /j,^^^ through Eq. (13b), i. e. 

f dk' 1 

P V (27r)3^/3(4°)-.(o))^^- ^^^^ 

At first order, from (13a) one has 

4" = 4"-9/(||4°V*"-«'). (15) 



with = l/exp[/3(eg^ - /i^^)) + 1]; from (13b) one then has 

dk' 1 



(27r)3 5(4V-/.«) 



+ 1 



(16) 



which fixes a new chemical potential The latter is then put back in (13a) 
together with e^ ^ in order to generate e? ^ and so on. 

Numerically the procedure is quite straightforward. However, it is of interest 
to analyze its diagrammatic content, since it helps to understand the different 
roles played by the various classes of diagrams at zero and finite temperature. 

Thus — following the steps outlined above for the solution of the HF equa- 
tions — the first order self-energy S^^^ is displayed in Fig. 2a: it embodies 
IJ,^^\ the chemical potential for the non- interacting system. together with 
fj,^^^ — the new chemical potential determined by the requirement of fixed 
density, Eq. (16) — is then inserted into the Dyson equation for the thermal 
propagator, namely 

g{k, Un) = g\k, UJn) + Q\k, UJn)^*{k, T)g{k, Un) (17) 
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Fig. 2. The thermal first order self-energy (panel a) and the perturbative expansion 
of the thermal first- iteration HF propagator (panel b). 



whose solution is (Fig. 2b) 



(18) 



Fig. 2b clearly illustrates that the chemical potential associated with the ex- 
ternal propagation lines or with those linking different self-energy insertions 
should now not be fj,^^\ but rather fj,^^\ in order to keep the density of the 
system constant. 

— * 

One then computes with the propagator (18) the new self-energy E^^(A;, T), 
which will embody the HF diagrams shown in Fig. 3a and will be used to fix 
a new chemical potential //(^). 

Then, solving the Dyson equation with {k, T) as a kernel will obviously 
lead to the new propagator shown in Fig. 3b. The procedure should be iterated 
until self-consistency is reached. Through this iterative procedure all the HF 
diagrams are generated and naturally organized in different classes. 



2. 1 The case of a zero-range interaction 



As a preliminary we address the HF problem at finite T considering the simple 
case of a zero-range interaction of pure exchange nature, namely we set in (9) 




(b) 



Fig. 3. Perturbative expansion of the first iteration irreducible HF thermal 
self-energy (panel a) and of the second-iteration HF thermal propagator (panel 
b). 



v{xi - X2) = Vi6{xi - X2), 
whose Fourier transform is of course just Vi. 

The HF equations (13) in this case reduce to: 

dk' 1 



(19) 



-9^1 
d^' 



(27r)3e'3(^fe'-'^) + 1 



(20) 



(27r)3e^(^E'-'^) + l' 



At zero-order = e^^"* and the chemical potential (a few terms of its 
expansion as a function of A;^ and T are quoted in Appendix A) is fixed by 
the requirement of yielding the correct density (see Eq. (14)). 



By inserting these results into Eq. (20) we get 



dA;^ 1 
k "-^ I (27r)3g/3(4°)-M(°)) + i 



^4°>-£v.p. (21) 
Again, the chemical potential has now to be redefined in order to keep the 
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(a) (b) (c) 

Fig. 4. Diagrammatic representation of the equation H*^^ = T,^^ 

density constant (see Eq. (16)). It is immediately seen that in the present case 
of a zero-range interaction the change of /i amounts to a mere constant shift, 
i.e. 

= /.(°) - lv,p. (22) 
Performing then a second iteration we obtain 



dk' 



f-mj 



(2t)%/3(4;'-^<") ^ J 

=4°'-^.p=4". (23) 



where KE})^ is here computed using the propagator (18), which, for a zero- 
range interaction, reads 

= ^ M°) QT/ M ^^^^ 

Since Eq. (23) clearly entails the identity 

/lE^i) = hE*j,^, (25) 

without any depence on T, we conclude that self-consistency is immediately 
achieved at any temperature for a zero-range force. In other words, for such 
an interaction the HE problem is trivial both at T = and at finite T. 

The diagrammatic content of Eq. (25), displayed in Eig. 4, helps in under- 
standing this conclusion: in fact, Eq. (25) holds valid because the diagrams on 

the right hand side of Fig. 4 cancel at any temperature order by order in the 
coupling constant Vi. Note that Vi appears both in the interaction lines and in 
the chemical potential jj,^^^: hence even the diagram (a) in Eig. 4 contributes 
to all orders in Vi. 
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To illustrate analytically these cancellations we evaluate diagram (a) to the 
order 0{Vi) using Eq. (22). We get 



-9 



-9 



Vi f dk 1 
Vi r dk 1 



-9/3yip/4 



J (27r)%;3(4"'-M('») ^ ^ 

^ , E^=i(-9/3yip/4)Vn! 



dk 



■4¥^+^'^^) (2.)3 



+ O (V,') . (26) 



From the above one sees that on the right hand side of Eq. (25) the term 
linear in Vi cancels with the left hand side, whereas the quadratic term cancels 
with the leading order term in Vi associated with the diagram (b). Indeed, 
computing the latter at the lowest order, that is with //(°) in all the propagation 
lines, one obtains 




Performing then the sums over the Matsubara frequencies with the usual con- 
tour integration techniques as in Eq. (7), one gets a contribution precisely 
canceling the term quadratic in Vi in Eq. (26). By the same arguments the 
cancellation is proved to hold at any order in Vi. In the following we shall 
refer to this occurrence equivalently as the cancellation theorem or the KLW 
theorem. 



2.2 The case of a finite-range interaction 



Wc now show that the cancellation mechanism previously illustrated occurs 
(under certain conditions of symmetry) also for a finite range potential, but, 
strictly speaking, in the hmit T — > only. 
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For this purpose let 

v{x,-X2)^V, — - — , (28) 

in \xi — X2\ 

which yields back the zero-range force in the A — > oo hmit and whose Fourier 
transform is 

(note that in the figures we shall express the range parameter A in MeV) . The 
first order self-energy is then easily obtained and reads 

V, A^ 1 r ,,^,^,x' + {k+hr 1 



S;,.(A;,T) = -9-^7 — TTT— / dA;iA;iln— . (30) 

\ / g _|_ 

In the T — > limit, where the Fermi distribution reduces to a ^ function, the 
integration in Eq. (30) can be done analytically, yielding 



E^,)(/c,r = o) 




kp — k A;/,' + A; 

arctan — ^ h arctan 



A 



In 



i + ikp-ky/x^ 



A 



(31) 



At second order one gets 



A^ + (A,, + A:,)2 e<'--"') 



1 



^ p + {k^ - k,y ^ g<^-.(°)) + ^^^^ 

The above can be obtained from Eq. (8), provided one inserts fj, — //^°^ in 
the innermost propagator and fj, — ji^^^ elsewhere, according to the previous 
discussion (see Fig. 3a). 

Eq. (32) contains a Fermi distribution and a factor which is proportional to 
the derivative of a Fermi distribution with respect to the energy e[?^ . These, in 

the hmit of T — > 0, yield a theta function, 6{^^^^ — e?^), and a 5 distribution, 

5(//*^^^ — e||^), respectively. Thus also E^2)) this limit, can be analytically 
expressed as 
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A 



A2+ i^^j2mfx(^yh+kF 
A2 + (^^2mi^W/h - kp 



(33) 



and is not vanishing. Nevertheless, also for finite range forces self-consistency 
is immediately achieved for T — > 0. 



The proof proceeds along the hnes foUowcd in the case of a zero-range inter- 
action. Indeed, since the chemical potential at first order is given by 



(0) ^..A^itir/ A 2kF ln(l + 4fc|/A2)\ 

we can compute, through the same power expansion illustrated in Eq. (26), 
the diagram (a) of Fig. 4 to the order O (V^), getting 



E^,^(^,r = o) 



^(1) ^ 



0) + 



9Vi 



1 2 



(27r) 



rnkf 
~2k 



X 



A 2kF ln(l + 44/A2) 

1 - — arctan — — H 

kp A Aklp/X^ 



i„m±m:+o(v»). 



\^ + {k- kpf 



(35) 



We thus see that the second term on the right hand side of the above equation 
exactly cancels the contribution of order arising from the diagram (b) of 
Fig. 4, whose analytic expression is given by Eq. (33) computed with = 
h'^kp/2m and replacing /j,^^^ with /j,^^^ in all places, i.e. in all the propagation 
hnes. This is the essence of the cancellation theorem. 



Somewhat surprisingly the validity of the latter actually extends over a re- 
markably large range of temperatures, as it appears from the numerical results 
reported in Subsection 2.3. 
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2.3 Numerical results 



We now present a few numerical results to assess quantitatively the role of 
the higher order perturbativc HF diagrams at finite T and the behavior of the 
chemical potential as a function of temperature and density. Wc also explore 
the sensitivity of our outcomes to the values of the parameters characterizing 
our schematic interaction. 

Specifically, we compute: 

1) the temperature behavior of the self-consistent HF chemical potential /xhf 
at a fixed density for various strengths and ranges of the force (Fig. 5); 

2) the size of the difference — between the successive iterations leading 
to the HF chemical potential, as a function of T to asses the impact of the 
higher order perturbativc terms on /xhf (Fig- 6); 

3) the temperature dependence of the first and second order contributions to 
the HF self-energy (Fig. 7). 

4) the temperature and density dependence of fx^^^ — //hf for a given interaction 
(Fig. 8); 

In Fig. 5 the HF chemical potential is seen to be simply shifted from the non- 
interacting value fj,^^^ by a constant amount over a wide range of temperature, 
as it happens for a zero-range interaction. 

Notably this shift turns out to be essentially given by the first order proper 
self-energy computed, for T = 0, at the Fermi surface. Since KE*^^ turns out 
to be remarkably stable with the temperature (see Fig. 7), wc infer that the 
cancellation theorem, strictly valid only in the T ^ limit, actually remains 
operative over a quite large range of temperatures, almost until the Fermi 
temperature. This is defined as 

keTp ^eF=^ (36) 

and one has ksTp = 34.48 MeV (/?f = 0.029 MeV"^) for the interaction in 
Eq. (28) with A = 140 MeV and Vi = 50 MeV fm''^ and ksTp = 26.7 MeV 
{Pf = 0.037 MeV"^) for the same interaction, but with Vi = 150 MeV fm^. 

Wc then conclude that for temperatures up to the proximity of Tp the chemical 
potential is 

AtHF = e^? + ^S^i)(A; = A;^., T = 0), (37) 

as follows by comparing the T — > limit of the self-energy shown in Fig. 7 
(heavy dashed fine) with the absolute value for large (3 (3.89 MeV) of the 
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Fig. 5. The difference — /xhf as a function of /3 for the cases of a force of 
zero-range (dotted) and finite range: A = 3000 MeV (dashed), A = 770 MeV 
(dot-dashed), A = 140 MeV (sohd). The hght lines refer to a strength 
Vi = 150 MeV fm^, the heavy ones to V-y = 50 MeV fm^. The density is p = 0.17 
nucleons/fm^. 




le-06 



P (MeV ) 



Fig. 6. The differences computed in HF at finite T, as function of (3 

in the case of a finite range interaction with A = 140 MeV and Vi = 150 MeV fm^. 
Note that the peak of the curves essentially corresponds to the Fermi temperature. 
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Fig. 7. nS^^) (dashed) and /iS^2) (solid) as a function of /3 with A = 140 MeV 
and Vi = 50 MeV fm^ for two different values of A:: k = (light lines) and 
k = kp = 1.36 fm~l (heavy lines). 




Fig. 8. fi^'^^ — /iHF as a function of /3 and p for a finite range force with A = 140 MeV 
and Vi = 150 MeV fm^. 
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heavy solid line in Fig. 5, which expresses /i^^^ — //hf- 



It is also evident in Fig. 5 that /xhf goes to /i^^^ for small values of /3 (classical 
limit), as it should, since here all the Feynman diagrams vanish. In fact, our 
exchange interaction is of purely quantum nature and, as such, is bound to 
vanish at large temperature. This behavior is also transparent in Fig. 7. 

The findings of Fig. 5 are complemented by those of Fig. 6, which conveys the 
information on the number of iterations required to achieve self-consistency. 
The figure shows that higher order iterations (i.e. diagrams) become signifi- 
cant for T approaching Tp, the temperature separating the classical from the 
quantum regime. This outcome illustrates the significance of Tp: in fact, a 
degenerate normal Fermi system, namely with a well-defined Fermi surface, 
lives in the temperature range T < Tp. Accordingly, the impact of the Fock's 
diagrams of order larger than one — that tend to disrupt the Fermi surface — 
become appreciable precisely at the Fermi temperature. Away from Tp their 
contribution is gradually disappearing both for T — > and T — > oo but, as 
previously discussed, for two radically different reasons, i. e. the cancellation 
theorem and the classical behavior (all the diagrams going to zero), respec- 
tively. Hence, for small T the first iteration provides an almost exact estimate 
of the self-energy, in spite of the fact that, as we can clearly see in Fig. 7, /iS*2), 
far from vanishing at small values of T, approaches the asymptotic limit given 



In concluding this Section we offer a global view of the behavior of /x'^^' — /ihf 
as a function of both f3 and p in Fig. 8, where the features of the chemical 
potential above discussed clearly appear. 



3 Propagation in presence of an induced symmetry breciking: the 
case of the external magnetic field 

In this Section we address the problem of the relationship between the Matsub- 
ara and Feynman theories at T = in the presence of a dynamically induced 
symmetry breaking. In this situation the two theories provide different results 
[5], the correct ones being given by Matsubara. We revisit this issue first for 
the well-known example of a non-interacting system of neutrons placed in an 
external constant and uniform magnetic field pointing into the 2;-direction, 
whose second-quantized hamiltonian reads 



by Eq. (33). 




2m 



] 



(38) 



where 




(39) 
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a and /3 being spin indices and jio the magnetic moment. 

Next we switch on the interaction among the neutrons: in this case, as we 
shall see, interesting and, to our knowledge, new aspects of the temperature 
HF theory emerge. 

In the interacting case the KLW theorem no longer holds for the static suscep- 
tibility (namely the one obtained with the single-particle propagator) . However 
the zero temperature polarization propagator leads to a dynamical suscepti- 
bility identical to the static one, obtained in the T — > Matsubara framework, 
in the limit of vanishing momenta. 



3.1 Static magnetization: the propagator for the non-interacting Fermi sys- 
tem 



We first briefiy consider the calculation of the free Fermi gas static magneti- 
zation M through the standard formulas 



at T = 0, and 



/3a 

{M) = ^1^0 J (lxY,<^PaQa0{xT,XT^ 
Pa 



(40) 



(41) 



at finite T, Q, being the large volume enclosing the Fermi gas. 



One could think to obtain the zero-tcmpcraturc propagator in Eq. (40) by 
solving the Dyson equation with the constant proper self-energy 



Kb = (-1) 



a+- 



i.e. viewing the external field as a perturbation (see Fig. 9). 
One would end up with the expression 



(42) 



Gap{k,Uj) = 5af3 



e{k - kp) ^ 9{kF - k) 



(43) 



which is diagonal in spin space, but not proportional to the unit matrix since 



^ka 



2m 



h 



(44) 



In Eq. (43) kp coincides with the Fermi momentum of the non- interacting 
system before switching on the magnetic field. 
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Fig. 9. The Dyson's propagator for a non-interacting Fermi system in a uniform, 
constant magnetic field. 



(a) 



(b) 



Fig. 10. The Fermi energy for a non-interacting Fermi system: in its ground state (a) 
and in the ground state in an external magnetic field (b). The two boxes represent 
particles with spin up and down, respectively. 



However, by inserting Eq. (43) into Eq. (40), it is immediately found that 



<M> = --;.o 



da; 
dk 



[ dk f 



dk 

da; 



duj 



e{k - kp) 



I3a 

eik-kp) 
e{kF - k) 



+ 



e{kF - k) 



Uk+ - IT] 



UJ — ujk- + iv Oi) — ujk- — iv 



0, 



(45) 



since the contributions of the two poles lying in the upper (Ima; > 0) complex 
energy plane cancel out [5]. 

This result is obviously wrong, since the propagator in Eq. (43) (diagrammat- 
ically displayed in Fig. (9)) does not correspond to an equilibrium state of 
the system (i.e. it does not correspond to a minimum of the energy): the true 
ground state of the system (see Fig. 10) is obviously unreachable perturba- 
tively, since the vertex in Eq. (39) cannot flip the spin of the constituents. 



The situation is different at finite T. Here the exact Matsubara propagator 
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reads 



lUJr, 



(46) 



and leads to the magnetization 



— * 



X 



/ 



(2^ 



-(0) 



(47) 



Clearly, Eq. (47) , together with the condition 



dfc 
(2^ 



+ 



(48) 



fixing the density of the system, in general can be computed only numerically. 
However, if the magnetic field contributes weakly to the energy eigenvalues, 
one gets the explicit expression 



(M) = -/.o / 



dfc 
(2^ 



n{ef + i,oB) - n{ef - i^oB) 
dk I dn 



,(0) 



-2ixlB 



k 



(27r)3 V de,- 

which, as it is well-known, in the T ^ limit, yields the magnetization 

Si^lpB 



(49) 



(M) 



and the susceptibility 



X 



2eF 

d{M) _ 3^1p 



(50) 



(51) 



dB 2eF 

namely the Pauli paramagnetism. Moreover, in the high temperature limit, 
the Curie law 



X 



(52) 



is attained. Thus, Matsubara theory, rather than perturbation theory at T = 0, 
yields the correct result. 
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The above findings relate to the structure of the thermal propagator, which 
sums over all the possible configurations weighted by the statistical operator 
and hence leads to the true, symmetry-breaking, ground state in the T — > 
limit. 

In fact, by computing the real time Green's function at finite T through a 
proper analytical continuation of the thermal propagator and then by taking 
the vanishing temperature limit, one gets 



Ga(3{k,Uj) = 6af: 



0{el - ep) e{eF 



e 



+ 



k' 



uj — ujka + ^ ^€f + ir] uj — ujka + fi- ^ep — irj 



(53) 



By comparing the above with (43), we see that an additional Fermi energy 
appears in the denominator and, most importantly, kp is replaced by the 
Fermi energy ep and A; by eS. Thus, the poles in the complex a;-plane of the 
propagators (53) and (43) (or the associated spectra) are different: hence the 
contributions they provide to the frequency integral no longer cancel out. 



3.2 Dynamic magnetization: the non-interacting linear response theory 



In this Subsection we briefly derive the magnetization of the free Fermi gas 
in the linear response theory at T = 0. It is well-known that this framework 
leads to the correct result for the Pauli paramagnetism, since it allows a weak 
external magnetic field B{x,t) to break the rotational symmetry of the non- 
interacting ground state even at T = 0, thus leading to magnetization. 

We start by recalling that turning on a perturbation H^^^ at time t = to, the 
fiuctuation of the vacuum expectation value of a generic operator O at time t 
reads 



6 < d{t) > = < A.t{t)\0\iJoAt) >-< ^olOlV^o > 

= £ dt' < V'ol [HTit'), OHit)] > . (54) 

In our simple example the perturbation is 

H'^\t) = -Ho J dxi^lix, t)a'i>H{x, t)B{x, t) (55) 

and the role of O is played by 

alj{x, t) = ip\j{x, t)a^ipH{x, t), (56) 
I ■00 > being the Fermi sphere. 
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Then, the magnetization along the z axis, namely the magnetic moment of 
the volume element, 

<M^{x,t)>= \XQ<a''{x,t)>, (57) 
is easily expressed as follows 

< M,(f , t)>= ^-iil j df dt'e{t - t') < 0o| [a|,(f , t), a|,(f', t')] 100 > B{x', t') 
^-l4 J dx'dt'U'''^{xt,x't')B{x',t'), (58) 

where the retarded spin-spin polarization propagator at zero-order 

im°J(a, x't') = e{t - t') < 41 [&Ij{x, t), alfix', t')] \4>o > (59) 

has been introduced. By Fourier-transforming (our system is translationally 
invariant), Eq. (58) becomes, for a; > 0, 

< Mz{k,uj) >= -nlTl^{k,uj)B{k,uj), (60) 

n°(A;,c<j) being the familiar free Fermi gas polarization propagator [6]. The 
above, in the c<j ^ 0^ and /c — limits, yields the magnetization induced by 
a static and uniform magnetic field. Since [6] 

hmRen»(i,0) = -j^ = --^, (61) 

one gets 

< >= --^^ilB, (62) 
thus recovering the correct susceptibility given by equation (51). 

In Eq. (61), the a; — > hmit, which imphes a static external field, should be 
taken before the A; — > one, which implies a uniform external field. 

3.3 Static magnetization: the propagator for the interacting Fermi system 

We now enlarge the previous analysis by switching on the neutron-neutron 
interaction 

^2 g-A|xi-S2| 

V{xi- X2) ^ ai- a2Vi —— — . (63) 

AtT \Xi — X2\ 

The first order self-energy splits now into two pieces that, at finite T, read 
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1 



+ 1 e + 1. 



(2^ 



Le 



+ 1 



/3(4°'+MoB-m) 



+ 1. 



(64a) 



and 



(27r)3 

dA?i 



i;(A; — ki) 



l3{ef'>-ixoB-ix) 

fci -1-1 



+ 



e '=1 
1 



e ''I 
1 



/3(6^°^-MoB-m) 

e 



The associated HF equations then become 



v{k — ki) 




+ 1 



1 



+ 1. 



+ 1 

(64b) 



9 



+ 1 



+ 1. 



+ 1 e 



+ 1. 

9 



+ 1 



e +1, 



(65) 



dfci 
(2^ 



,e '^i +1 



e ^=1 + 1, 



dfci 
(2^ 



,e + 1 



+ 



6*1 + ly 



From Eqs. (64) and (65) it appears that, owing to the presence of the magnetic 
field the self-energy not only is no longer proportional to the unit matrix 
in spin space, but it acquires, beyond the exchange contribution, a direct one 
as well, an occurrence with far reaching consequences. 

We have solved the system (65) numerically, searching for the magnetization in 
Eq. (47) and the HF chemical potential //hf- Our results for //hf are displayed 
in Figs. 11 and 12. For purpose of illustration wc have chosen the very large 
value of 10^^ tesla for B (on the surface of neutron stars B is estimated to 
range from 10^ to 10^'' tesla [7], whereas in the interior its upper limit is set 
at a few times 10^^ tesla 



Remarkably, in the presence of a magnetic field the contributions to fi^p arising 
from the higher iterations (Fig. lib) are seen to stay quite constant with the 
temperature till very small T, unlike the case with no symmetry breaking 
(Fig. 11a), where they reach a maximum for temperatures close to Tp. 
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Fig. 11. The difference between the successive HF iterations /x^*-* — fi 
function of /3 in the case of a finite range interaction with A = 140 MeV and 
Vi = 50 MeV fm^ in absence of magnetic field (a) and with a background magnetic 
field B = 10^^ tesla (b); p = 0.17 fm~^. In panel (b) the modulus of the differences is 
displayed: note the alternation in sign that implies an oscillating approach to /xhf- 
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Fig. 12. As in Fig. 11, but with Vi = — 50 MeV fm^. Note that the differences are 
now constant in sign, except the first one: hence the approach to /Uhf is smooth. 

In the proximity of T = 0, the higher order diagrams contributing to /xhf, 
rather than cancehng each other, now contribute more than they do at high 
T. This is a beautiful manifestation of the failure of the cancellation theorem 
in the presence of a symmetry breaking. 



Of significance is also the alternating behavior of the iterations for an antifer- 
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romagnetic interaction: those of odd order have a positive sign, those of even 
order a negative sign (note that in Fig. lib we display the absolute value of 
As a consequence the approach to the HF self- consistent solution 
is no longer smooth, as in the symmetric case, but now the successive itera- 
tions oscillate around the HF solution until they stabihze at the latter value 
when self-consistency is reached. 

In contrast, it is remarkable that for the ferromagnetic coupling the sign of the 
successive iterations is constant (positive), except for the first one, which is 
negative. Hence, in this case the self-consistent solution is smoothly reached. 
The situation is reminiscent of the physics of an hydrogen atom in an external 
magnetic field, the competition between the Coulomb and the magnetic field 
now being replaced by the one between the nucleon-nucleon interaction and 
the magnetic field. If this opposes the action of the nucleon-nucleon forces, 
then it makes it harder to reach self-consistency. 



3.4 Dynamic magnetization: the interacting linear response theory 



As previously seen for the case of a non-interacting Fermi system, also when 
the interaction is switched on we can obtain the system's magnetization and 
susceptibility in the framework of the linear response theory at T = using the 
polarization propagator, a tool set up with the two-fermion Green's function. 

We start by computing the spin-spin retarded polarization propagator needed 
in the HF theory. In the present case, owing to the nature of the interaction 
in Eq. (63), only the Fock term is non- vanishing, so we shall actually compute 
H^. At T = its real part reads 



ReU^{q,u;) = -V J ^^(| k + q\ -kFMk, - k)-^-^^, (66) 



where 

^kq = ^Uq + ^(1) i^' + ^-'^k- ^(1) (^)- (67) 

The above expression cannot be computed analytically. Yet, in the case of the 
interaction in Eq. (63), a quite accurate analytic approximation to Eq. (66) 
can be obtained by expanding the self-energy in Eq. (31) (divided by 3, since 
we dropped the isospin degree of freedom) up to and including terms in k'^; 
one gets 

S^^(fc) = A + (68) 

with 

^ = ^Lctan^-^) (69) 
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and 



B 



(70) 



— * — * 

Then, computing S*(A; + q) — S*(/c) in this approximation and inserting this 
result into Eq. (66), we get for RcIT^ an expression identical to the free one, 
but for the replacement of the nucleon mass m with an effective mass m* given 
by 



m 
m 



1 + 



h^7r^{kj, + x^y 



(71) 



Next, letting as before first a; ^ and then g — > 0, we obtain for the magne- 
tization and the magnetic susceptibility the expressions 



Z 6rr 



(72a) 



and 



X 



(72b) 



respectively, where the Fermi energy now reads €*p — h'^kp/2m*. Since, for the 
interaction in Eq. (63) with Vi > 0, one has m* < m, we conclude that for an 
antiferromagnetic force the HE mean field lowers the free magnetization (or 
the susceptibility). 



However it is of significance that Eq. (72a) does not coincide with the T 
limit of the magnetization 







(M) = /io 



dk 



e^^'l-^^ + 1 



e^^^'t-^^ + 1 



(73) 



ei", and ^ being the HE single-particle energies and chemical potential 
obtained by solving Eq. (65). To achieve the accord between Eqs. (72a) and 
(73) it is necessary to go beyond HE. This will be done in the next Section. 



Here we display in Eig. 13 the temperature behavior of the relative magneti- 
zation M (yielding the fraction of spins anti- aligned to the magnetic field, i.e. 
M =< > /lAiolP = {N^ - N'^)/{N^ + N'') with obvious meaning of the 
symbols) associated with Eq. (73). Note the decreasing of M as T increases 
and the recovery at large T of the Curie value, the interaction among neutrons 
becoming irrelevant at high T . Moreover the values of M at T = are lower 
(larger) than the free one for an antiferromagnetic (ferromagnetic) force: they 
will be shown to coincide with the predictions of the linear response theory in 
a frame extending the HE one. 
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Fig. 13. The HF magnetization in an external magnetic field as a function of 
temperature for B = 10^^ tesla and for different values of Vi; A = 140 MeV and 
/> = 0.17 fm-3. 

4 Beyond HF 



In the framework of the linear response theory the natural extension of the 
HF approach should embody the RPA (random phase approximation) corre- 
lations. To start with, we neglect the exchange terms, i.e. we make use of the 
simpler ring approximation. 

This yields for the spin-spin polarization propagator of a homogeneous infinite 
system the expression 

ff-(g,a;) = , (74) 



from where the magnetization 



<M>=- lim lim Ren"'^s(g' a;)^2^ = ^^(^5^^ 2^ 

l + yi(3p/2e^)^° ' ^ ^ 

which is linear in B, is deduced. Thus, for the magnetic susceptibility in ring 
approximation, one gets 

ring _ 3p/2e;^ 2 

showing that the ring correlations, for an anti-ferromagnetic (Vi > 0) nucleon- 
nucleon interaction, lower the susceptibility of the free Fermi gas, thus acting 
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in the same direction of the HF mean field. 



The quenching of the free Fermi gas susceptibihty x^^^*^ ^i^^g approximation 
relates to the existence of a spin collective mode, called a magnon, which 
makes harder the excitation of the system. As it is well-known, the energy of 
the magnon is found by solving the equation 

1 - Vi^^ReU'^iq^u) = 0, (77) 
g"^ + A"^ 

that — in the limit a; — >• 0, g — >• 0, but keeping now the ratio x = {u}/q)m* / hkp 
constant and greater than one — is fulfilled by the phonon-like dispersion 
relation uj — esq, the zero-sound velocity cs being fixed by the equation [6] 

with 

= -xln^^-l. (79) 

To assess the role of anti-symmetrization one resorts to the full RPA, which 
includes the exchange diagrams beyond the ring ones. This we do with the 
method of the continued fractions [9,10,11], which, when truncated at the first 
fraction, yields for the RPA polarization propagator the following expression 

n-^(.» = p ""-^^^ ' (80) 

1 - V,^^U^q,u;) - n^P{q,u;)/U^q,u;) 



q' + X' 



where 



nr«(^,^) = -/ (^^A'^) / ^,E{k,,q)v{h-k,), (81) 



being 



E{k, q) = ^ 



9{\k + q\- kF)9{kF - k) 9{kF -\k + q])9{k - kp) 



, (82) 



is the exchange diagram associated to the first order ring contribution. The 
expression in Eq. (80), while exact for a zero-range force, remains a good 
approximation for typical nuclear finite range interactions too. 

One can verify that the magnetization deduced from Eq. (80) embodies a 
susceptibility very close to the exact RPA one, since the RPA problem can be 
analytically solved for an infinite system in the limit of vanishing frequencies 
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and momenta, where it merges with the Landau's quasi-particle theory [12,13]. 
One gets (see Appendix B) 



RPA ^ 'ip/'^^t 2 coo\ 

^ l + yi(3p/2e^) + |yi(3p/2e:^)(A/2M2 1n[l + (2A;^/A)2]^°' ^ ' 

which, for a zero-range force, becomes 

RPA ^ ?>pl2eF , . 

Xo-range iJ^Wip/Aep ^ ' 

an expression identical to the one deduced from (80). Eq. (83) shows that, in 
the anti-ferromagnetic case, the anti-symmetrization reinforces (mildly) the 
quenching of x^''^^ induced by the ring correlations. 

Note that in RPA the impact of the exchange diagrams decreases as the range 
of the force increases, as it should: this result follows by comparing Eq. (83) 
with Eq. (84). 

Wc arc now in a position to find out the polarization propagator yielding (as 
long as the response of the system is linear) the HP Matsubara magnetization. 
For the sake of simplicity we start by considering a zero-range antiferromag- 
nctic force of strength Vi = 50 MeV. In Pigs. 14 and 15 are displayed versus 
B various relative magnetizations as obtained from the Matsubara propagator 
and from the zcro-tcmpcrature linear response theory. In Pig. 14 wc employ a 
zero-range force, hence the HP mean field does not contribute to the polariza- 
tion propagators; in Pig. 15 a finite range force. The figures clearly show that 
the magnetization obtained from n^^'^ coincides with the one given by the 
HP thermal propagator Gffp for T ^ until ~ 3 (in units of 10^^ tesla), 
which sets the limit of validity of the linear response framework. Por B larger 
than B the magnetization increases much less rapidly than linearly until it 
saturates when the system is fully magnetized, an expected behavior correctly 
reproduced by the thermal HP theory. We thus conclude that in order to 
obtain the magnetization and, in general, any mean value using the single- 
particle thermal propagator Q^p, it is necessary to sum up an infinite series 
of loop diagrams in the particle-hole polarization propagator, namely those of 
RPA. This correspondence extends also to the pure ring and ladder polariza- 
tion propagator as illustrated in Table 1 (with n''^'^^'^'' we mean the propagator 
summing up only the particle-hole exchange diagram). 

Also worth noticing is that the range of validity of the linear response theory 
appears to be only moderately affected by the many body scheme employed. 

Turning now to comment on the finite range force, we observe that the slopes of 
the curves in Pig. (15) are somewhat reduced with respect to those in Pig. (14). 
Indeed, in the former the fermionic lines in the polarization propagator are 
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Fig. 14. The magnetization as a function of B for Vi = 50 MeV fm^ and zero range 
atr = {p = 0.17 fm~^). The dotted hnes display the results in the linear response 

framework (using 11), the soHd Hnes in the thermal self-consistent calculations at 
r = (using Q^). Starting from above the pairs of curves correspond, respectively, 
to the following cases: W^/G^, U^'"^'^"' /Q^ , n^'^s/^H' ^^^^/Gip- 




Fig. 15. The magnetization as a function of B for Vi = 50 MeV fm^ and 
A = 140 MeV at T = (/9 = 0.17 fm^^). The dotted hue displays the result in 
the free linear response framework (II'^), the upper dashed line in the linear re- 
sponse framework with Fock correlations (XIp) and the lower dashed line with Fock 
and RPA correlations (11^^^); the upper solid line represents the results of the ther- 
mal self-consistent calculations at T = including only Fock correlations {Qp), the 
lower solid line includes HF correlations (Q^p)- 
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Table 1 

The zero-temperature polarization propagator required to obtain the same mag- 
netic susceptibility deduced (in the T — > limit) from the Matsubara free (Gq), 
Hartree (G^), Fock (^p ) and HF (G^p) thermal single-particic propagator. Clearly 
the many body diagrams summed up to infinite order to get the appropriate 11 are 
the ring, the ladder and the RPA ones, respectively. 

dressed by a dependent Fock self-energy and this, as previously shown, 
reduces the susceptibility for an antiferromagnetic interaction. Furthermore, 
as expected, Hp™^ and FIp^^ tend to become very similar as the range of the 
force increases. 

Thus we have reached the result that the statistical average required in com- 
puting is performed over the states entering into the spectral represen- 
tation of the RPA polarization propagator with fermionic lines dressed by a 
Fock self-energy. 

In concluding this Section we display in Fig. 16 the frequency behavior of the 
magnons, both in ring and in RPA, with and without the HF mean field, for 
an interaction with Vi = 480 MeV fm^ and A = 140 MeV. The hardening of 
the magnon mode stemming from the HF mean field and (to less extent) from 
anti-symmetrizing the ring diagrams is clearly apparent in the figure. 

We display also in Fig. 17 the system's spin response functions (proportional 
to the imaginary part of H(^, a;)) both in ring and in RPA and both for an 
antiferromagnetic and a ferromagnetic coupling at g = 1 MeV/c. For sake 
of illustrating the impact of the pure ring and RPA correlations on the re- 
sponse we ignore in the figure the action of the Fock mean field (namely, we 
compute Eqs. (74) and (80) replacing H^ with H°, or, in other words, we set 
m* — m). The values chosen for the coupling are Vi — ±150,480, —239 and 
—223 MeV fm^, whereas the range parameter is A = 140 MeV. In the figure 
it is clearly seen that: 

a) for an anti-ferromagnetic coupling, as expected, the magnons are standing 
out above the particle-hole (ph) continuum, the more so the larger Vi is. 
Correspondingly, the depletion of the ph continuum increases with Vi. Note 
also that the impact of anti-symmetrization is modest, as expected owing 
to the long range of the force: the more so, the larger Vi is; 

b) for the ferromagnetic coupling the response is of course enhanced, the en- 
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Fig. 16. The dispersion relation for the collective modes (magnons) in ring (dashed) 
and RPA (solid), with (upper lines) and without (lower lines) the HF mean field 
for an antiferromagnetic interaction with Vi = 480 MeV fm^ and A = 140 MeV 
{p = 0.17 fm~^); the dotted line represents the upper border of the ph response 
region, q'^/2m + qkF/m. 
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Fig. 17. Spin response functions at g = 1 MeV/c, p = 0.17 fm~^ and for a fi- 
nite range interaction with A = 140 MeV. The dotted lines represent the free re- 
sponse, the dashed and solid lines the ring and RPA ones, respectively. In panel 

(a) we employ antiferromagnetic couplings: Vi = 150 MeV fm^ (upper curves) and 
Vi = 480 McV fm^ lower curves; the spikes represent the ring and RPA collective 
states: the pair at higher energy corresponds to Vi = 480 MeV fm^, the other one to 
Vi = 150 MeV fm^ (note that in the latter case the ring and RPA excitations energies 
are very close). In panel (b) we employ ferromagnetic couplings: Vi = —150 MeV fm^ 
(lower curves), Vi = Vl'"'^^ ^ -239 MeV fm^ and Vi = V^^^^ ^ -223 MeV fm^ 

(note that the ring response at Vi = V^"crft ^^'^ RPA one at Vi = V^^^ practi- 
cally overlap each other). For the sake of simplicity, the HF mean field has not been 
included. 
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hancement becoming however dramatic for Vi = —239 MeV fm^ in ring 
and for Vi = —223 MeV fm^ in RPA. These values of the couphng are 
those yielding a divergent ring and RPA susceptibility, respectively, when 
m* — m. We shall return on this point in the next Section. 



5 The spontaneous symmetry breaking 



From the explicit formulas of x"^^"^ and x^^^ obtained in the linear response 
scheme one sees that both susceptibilities diverge in the case of a ferromagnetic 
coupling in correspondence to the following critical values of the strength of 
the force: 

..ring _ 2e^ ring _ ^ ^ + ^ ) /or„\ 

l/l,crit - 3p ^ ^l>crit - + + 1 ^^^^^ 

and 



2et 



3p 1 + |(A/2A;,.)2 In [1 + (2A;,./A)2] 

'^l,crit — \4 , 0\9 , 1 , , \9.\9J\ /r,\9 1„[-t , /'0/\\91 \^^'^) 



mkp A4 + 3A2 + 1 + 1(1 + A2)2(A/2)2 ln[l + (2/A)2] 

with A = X/kp- The above expressions have been deduced making use of 
Eq. (71), which yields the neutron effective mass m* in terms of Vi, A and kp- 

The behavior of Eqs. (85a) and (85b) versus \ = X/kp is displayed in Fig. 18: 
note that the two curves coincide for an infinite range force (A = 0) whereas 
their difference is the largest for a contact interaction (A = oo). For A = 140 
MeV and kp = 1.71 fm-^ one has V^t = -239 MeV fm^ and V^^^^ ^ 
—223 MeV fm^: these are the values previously employed (see Fig. 17). 

In this Section, for sake of simplicity, we shall consider a zero-range force 
(A = oo). In this instance, always with kp = 1.71 fm~^, which corresponds 
to a neutron density p — 0.17 fm~^, one gets V\"crft — —239 MeV fm^ and 
^ifcru = — 159 MeV fm^, respectively. We remark that the value of V™^ for a 
zero range force coincides with the one of Ref. [3], which in turn is identical 
to what one gets in the Hartree theory, as we shall later see. 

The values of the couphngs in Eq. (85) signal the occurrence of a spontaneous 
symmetry breaking: correspondingly, a phase transition takes place in the 
system from a phase fully symmetric to a phase where the rotational invariance 
is broken in spin space. Hence, the vacuum appropriate to this new phase 
should be characterized by two Fermi momenta kp and kp, associated with 
neutrons with spin up and down, respectively. 
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Fig. 18. The behavior of the critical interactions V^i™ft (dashed) and V^^^ (solid) 
versus A = X/kp in units of — Tr^^^/mfc^. The arrow corresponds to A = 140 MeV 
and kp = 1.71 fm~^. 

In this situation the free Fermi propagator, while still diagonal in spin space, 
is no longer proportional to the unit matrix, but becomes 



Gll{k,u) 



\ 



G''±{k,uj) J 



(86) 



where 



G°++{k,uj) 
G'^ik.uj) 



9{k-k+) ^ 9{kt-k) 



eik-kp) ^ e{kp-k) 



and ujj^ = hk /2m — 



u — oj^ + ir) u — oj^ — ir] 

(0) 



(87a) 
(87b) 



With the above propagator the first order momentum independent self-energy 
is easily computed and, as expected, splits into two terms, namely 



HE 
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(ji'p kp ^ ~l~ 2kp ^ 

(jip kp ^ (j^p ~l~ ^kp ^ 



(88a) 
(88b) 



of obvious physical meaning. 
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In Eq. (88) both the direct (first term on the right hand side) and the exchange 
(second term) contributions are neatly separated out. Note that the direct 

term is the same in both and S , but for the sign, and moreover it 

vanishes, as it should, for kp = kp. 

Since for a homogeneous infinite system the HF problem is trivial at T = 0, 
even in a broken vacuum, then the exact HF fermion propagator is easily 
obtained by replacing Eq. (87) with 

e{k - kt) ^ e{kt - k) 

Oi) — (uj^ — SVikp^ /Qn'^h) +ir] u — (u^ — 3Vikp^ /dn^H) — irj 

(89a) 

e{k - kp) ^ e{kp - k) 

Oi) — (a;g — SVikp^ /dn'^h) +ir] uj — (u^ — SVikp^/dn^h) — irj 

(89b) 

The above, to be viewed as a generalized HF fermion propagator in a vacuum 
not invariant for spin rotations, embodies the HF single-particle energies 

and allows to compute the HF energy per particle in terms of the relative 
magnetization M = {N^ - N^)/{N^ + N^) and of the density p = {N^ + N^)/n 
according to 



G+V (fc,a;) = 

G ' [k,U!) ~ 



E 
N 



dk 
(2^ 



,(0) 



2 '67r2 



e{k+ - k) 



+- 



,(0) 
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Ak 



(27r)3 
1 + Mf'^ 



(0) 3 kp?> 
^ 2 '67r2 



e{kp - k) 



(1 - Mf'^] - '^Vip{l - M^) 



(91) 



Here N"^ /VL = kp'^/Gir'^ and N^/il = kp^/Qn"^ correspond to the spin up and 
spin down densities, respectively. Likewise, one can easily write down the 
expressions for ei and E/N m Hartree and in Fock approximation separately. 
We remark that Eq. (90) coincides with the single-particle energies obtained 
by solving Eq. (65), with a zero-range force, in the hmit T — > 0, S — > 0. 

Now, for the system to be in equilibrium the two chemical potentials (Fermi 
energies at T = 0) of the spin up and spin down fermions should be equal. 
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Hence, in HF the equation 

h'kf _3V, s_h'k-/ _3V, s 

should be fulfilled and similar equations should hold in the Hartree and in the 
Fock approximation as well. Notably all the three cases can be compactified 
into the single equation 

(1 + M)2/3 - (1 - M)2/3 = -a^^M, (93) 

eF[0) 

with 

a = 3 in HF 

a = 2 in Hartree (94) 
a = 1 in Fock. 

In the context of solid state physics Eq. (93) is usually referred to as Stoner 
equation. Of course, the values in Eq. (94) for a are valid only for a zero-range 
force. 

It is a remarkable occurrence that one arrives exactly at the same Eq. (93) by 
minimizing the energy per particle in Eq. (91) with respect to the magneti- 
zation. This outcome is a consequence of the Hugenholtz-Van Hove theorem 
[14] , which remains true even in a broken vacuum. 

Now, the Stoner equation in Eq. (93) admits real solutions only for M in the 
range < M < 1. For a given density p, which also fixes e^'', the case M = 1, 
corresponding to a fully magnetized system, leads to the critical value for the 
coupling 

02/3 (0) 

= — ^. (95) 

a p 

This expression cannot be derived in the linear response framework and with 
a = 2 (Hartree approximation) is identical to the mean field value quoted in 
Ref. [3]. 

Instead, in the case M — 0'^ (incipient ferromagnetism) Eq. (93) leads to 

4 e^°^ 
6a p 

This result, when a — 2 (Hartree approximation), coincides with the linear 
response value in ring approximation in Eq. (85a), whereas, when a — 3 (HF 
approximation), coincides with the linear response value in RPA in Eq. (85b). 

It is clear that with the propagator in Eq. (89) it would now be possible 
to study the system's response to an external probe when the vacuum is 
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broken and the associated collective modes. These split into excitations along 
the z-axis (assumed as the direction along which the spontaneous symmetry 
breaking has taken place) and along the directions orthogonal to it (Goldstone 
modes). This topic is presently addressed in a work in progress. 

We anticipate, however, that we expect to get in this case a collective mode 
characterized by a quadratic — not linear as in Eq. (77) — dispersion relation, 
which is the one actually exhibited by ferromagnetic crystals. Our guess is 
supported by the findings of Ref. [15], where, for a system (asymmetric nuclear 
matter) displaying a spontaneous symmetry breaking in the isospin space, a 
collective mode with a quadratic dispersion relation has indeed been found. 



6 Conclusions 

In this work we have reexamined the issue of the connection between the 
Matsubara theory at finite T and the Feynman one at T = addressing a 
specific example, namely an infinite homogeneous system of nucleons, and a 
specific theoretical many-body framework, namely the HF mean field. 

Our focus has been on how the theorem of Kohn, Luttinger and Ward is 
actually implemented. In fact, the KLW theorem assures the identity between 
the T = Feynman theory and the T ^ limit of the Matsubara one (at 
least for systems made up with spin one-half fermions) when no symmetry 
breaking occurs. 

Yet, one would like to know the actual temperature behavior of the HF di- 
agrams of order higher than one: these in fact rigorously vanish at T = 
— for a static interaction — when computed as Feynman diagrams, but not 
when computed as T — > Matsubara diagrams. Thus the theorem is realized 
through a cancellation, which becomes "complete" at T = and "minimal" 
at the Fermi temperature Tp. Hence, here is where the HF problem for an in- 
finite system becomes nearly as complex as in finite systems: indeed, reaching 
self-consistency at Tp requires to account for many diagrams beyond the first 
order ones or, equivalently, to perform many iterations when searching for a 
solution of the non-linear, non-local HF equations. Just how many it depends, 
of course, on the density and on the interaction and, concerning the latter, 
more on its range than on its strength: in fact, for a contact force the HF 
problem in an infinite system at any T remains as trivial as it is at T = 0, 
namely it simply amounts to first order perturbation theory. 

It becomes trivial as well at large T, at least when the interaction is of pure 
exchange (and hence of pure quantum) nature, which is the one we confine 
ourselves to consider, for sake of illustration, in this paper: large T indeed 
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means classical physics. In addition, the impact of any interaction is no longer 
felt at large T, where the kinetic energy dominates. 

When a symmetry is broken, so it is the KLW theorem. In the second part of 
this work we have addressed this issue again, aiming at exploring more closely 
how this occurrence is emerging in a specific case, meant to be illustrative 
of the phenomenon in general. We have thus studied a system of neutrons 
placed in an external magnetic field B, which dynamically breaks the rota- 
tional invariance of the system's Hamiltonian in spin space. For the purpose 
of clarifying the physics we choose the strength of B very large. 

Three basic features emerge from our analysis. First, although our two-body 
potential is of purely exchange nature — thus entailing only a Fock term in 
the interaction matrix elements — the presence of B restores the presence of 
a Hartree term as well, whose relevance dependes upon the observable one is 
dealing with. Thus, in our example, the role of the induced Hartree term is 
minor as far as the chemical potential in concerned, but substantial in the 
magnetization. 

Secondly, as T — > the successive iterations leading to the self-consistent field 
no longer cancel each other — as they would were the KLW theorem valid — 
rather they yield a net contribution, actually as large as it is at Tp, which stays 
pretty constant in the temperature range from TptoT — 0. We have verified 
the above outcome for the case of the chemical potential in the presence of an 
unrealistically large magnetic field; different values of B would affect the size 
of the contribution, but not its existence. 

Finally, it turns out that the approach to the self-consistent solution is much 
dependent upon the sign of the exchange interaction: indeed, we have found 
that a ferromagnetic coupling among the neutrons swiftly and smoothly leads 
to the HF mean field, whereas an antiferromagnetic coupling not only requires 
more iterations before reaching self-consistency, but also entails an approach 
to the latter which "oscillates" around the self-consistent solution. 

Given that in the presence of a dynamical symmetry breaking the Feynman 
and Matsubara theories in the T — > limit provide different results and that 
those of Matsubara are the correct ones, the question naturally arises: is it 
possible to obtain the right expectation values for the observables also in the 
framework of the zero-temperature theory? The answer to this question is 
positive, provided the appropriate propagator is used at T = to compute 
the observables. 

Thus, in the example we have treated, the chemical potential and the magneti- 
zation obtained from the HF single-particle thermal propagator in a magnetic 
field 

^HF can be obtained as well in the framework of Feynman theory at T = 
using the polarization propagator n(g, o;). The latter should however embody 
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the resummation of an infinite number of loops: specifically, the results ob- 
tained from Q^p in the T limit are recoverd in the T = framework 
using n^^"^. Actually this correspondence applies to any infinite subset of 
many-body diagrams, in the sense, for instance, that corresponds to 
and n'^^^"^ (the steps of the ladder being between a particle and a hole) corre- 
sponds to Qp. Of course, in the absence of any interaction the correspondence 
is between 11° and Qq. 

However, the above statements are vahd only in a limited range of values of 
B, namely as long as the linear response theory is tenable. It is worth noticing 
that we can actTially compute this range by comparing the predictions (c. g., 
on the magnetization) of the single-particle thermal propagator in the T ^ 
limit with those of the polarization propagator. It turns out that the range of 
applicability of the linear response theory is not much affected by the many- 
body scheme employed. 

A remarkable feature of the linear response theory lies in its ability to predict 
the onset of a phase transition in the system, in our case of the ferromagnetic 
phase, as a function not of the temperature, but rather of the strength of the 
coupling constant Vi viewed as a control parameter entering into the system's 
Hamiltonian. In fact, expressing the magnetic susceptibility of the system via 
the real part of 11 (g*, uj) and searching for the poles of the latter in the variable 
Vi, one finds [17], both in RPA and in ring approximation, the critical values 
^icrit ^i?crft whcrc the phase transition starts to occur if the coupling is 
ferromagnetic (incipient ferromagnetism) . 

Indeed, as it is well known, the spontaneous symmetry breaking is heralded 
by a dramatic increase of the isothermal susceptibility and by a spectacular 
enhancement of the spin response function at wavevectors small enough to 
encompass the long range spin order taking place in the system. Note that 
^i^rit ^"ctH substantially different for a contact force (by a factor 2/3, 
^/crft being the larger one, showing that the phase transition is made easier 
by antisymmetrization), while they become closer and closer as the range of 
the interaction increases, as it should be. 

For an antiferromagnetic coupling the above does not occur because, as already 
anticipated, the Fermi gas, even in the absence of the interaction, lives in the 
antiferromagnetic phase. Hence, it should display Goldstone modes with a 
linear relationship between frequency and wavevector, which indeed we have 
found. 

Naturally, one would like as well to explore the Goldstone modes in the fer- 
romagnetic case: this we have not done, but we have paved the way to their 
investigation using the technique of the anomalous propagator. In fact, in the 
ferromagnetic phase the system develops two Fermi momenta, kp and kp, and. 
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as a consequence, the single-particle propagator, while still diagonal in the spin 
indices, is no longer proportional to the unit matrix, but splits into the two 
components G^^'^ and G^^^. These we have computed and, with their help, 
we have computed as well the system's energy per particle in terms of the mag- 
netization M. By minimizing the energy with respect to M or, equivalently, 
by requiring identical chemical potential for the two spin species of particles, 
we have found the critical value of the coupling constant, V^^cri^) yielding a 
fully magnetized system, a result impossible to reach in the standard linear 
response theory framework. Worth mentioning is that the HF value of V^^^^'^ 
is lower than the corresponding values obtained in the simpler schemes of the 
Hartree's or Fock's theories. 

The above results clearly would not be attainable in the scheme of the Feyn- 
man perturbative theory at T = 0, even using the polarization propagator, 
which of course illustrates that in the presence of a spontaneous symmetry 
breaking the KLW theorem does not hold. 

In conclusion, it appears worth remarking that the free Fermi gas, a system 
to be viewed at T = as antiferromagnetic owing to its spin-zero wave func- 
tion, when the interaction is switched on displays both in the ferromagnetic 
and in the antiferromagnetic phases the same magnons found in crystals, in 
spite of being a continuous system. This recognition offers the opportunity to 
estabhsh a contact between the computations performed on the lattice and 
within advanced many-body frameworks based on the Fermi gas model, for 
example in connection with the Heisenberg Hamiltonian. Not however with 
the Ising model, which, although dealing with a physics not so far from the 
one dealt with in this paper, actually does not display the variety of collective 
(Goldstone) modes showing up in the present system. This occurrence reflects 
the fact that what is broken is a discrete symmetry in the Ising model, a 
continuous one in our infinite, homogeneous system of neutrons. 
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In this Appendix we recall the power expansion for the chemical potential of 
a free Fermi gas as a function of kBT/ep to the order 0{kBT /epY ■ 

Starting from the Sommerfeld's expansion of the density [16], 
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one solves the equation 
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iteratively. At third order in (the zero order being /j, — ep) one obtains 
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, (A.3) 



the last term on the right hand side not being quoted in the literature to our 
knowledge. 
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As shown in Rcf. [12], the RPA polarization propagator — approximately 
expressed by Eq. (80) in the first order of the continued fraction expansion — 
is exactly given in the long wavelength, low frequency limit by the formula 



lim n(a;) 



m*kr 



(277/1)2 1 + Mo' 



(B.l) 



with X — {uj/q)m* /hkp. In Eq. (B.l), uq is the amplitude of the I — partial 
wave entering into the expansion of the ph matrix element 



(B.2) 



which is vahd close to the Fermi surface (in general Wph(C) CO J^^Y have an 
arbitrarily large number M of components). 

In Eq. (B.2) ^ = k-q and ^' = k'-q, k and k' being the directions of the momenta 
of the initial and final hole, respectively, and q the direction of the momentum 
transfer. Moreover, in Eq. (B.2) the integration is over the azimuthal angles 
of k and k' and the Pi are the Legendre polynomials. 

Now, with the simple interaction we use, the a — 1 (spin one) ph matrix 
element reads 



V^^{k,k',q)^2V, 



A 



+ Vi 



X^ + \k-k'\'^' 



(B.3) 
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where the direct and the exchange contributions are exphcitly displayed. The 
direct one, easy to deal with, leads to the ring approximation expression in 
Eq. (76). Concerning the exchange term, from Eq. (B.2) one finds 



2^p^(e,0 = Vi(l + 4 



A4 



+ 4- 



A2 



again with A = X/kp- 

By comparing then Eqs. (B.2) and (B.4) one finally gets 
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{p = kp/3n^ being the density of our homogeneous neutron system), from 
which Eq. (83) is obtained. 
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